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Abstract. 12 CO J = 1-0, 2-1, 4-3, 7-6, and 13 CO 1- 
0, 2-1, and 3-2 line emission was mapped with angu- 
lar resolutions of 13" - 22" toward the nuclear region 
of the archetypical starburst galaxy M82. There are two 
hotspots on either side of the dynamical center, with the 
south-western lobe being slightly more prominent. Lobe 
spacings are not identical for all transitions: For the sub- 
millimeter CO lines, the spacing is ~15"; for the millime- 
ter lines (CO J = 2-1 and 1-0) the spacing is ^26", indi- 
cating the presence of a 'low' and a 'high' CO excitation 
component. 

A Large Velocity Gradient (LVG) excitation analysis 
of the submillimeter lines leads to inconsistencies, since 
area and volume filling factors are almost the same, re- 
sulting in cloud sizes along the lines-of-sight that match 
the entire size of the M82 starburst region. Nevertheless, 
LVG column densities agree with estimates derived from 
the dust emission in the far infrared and at submillimeter 
wavelengths. 22" beam averaged total column densities 
are 7V(CO) - 5 10 18 and N(R 2 ) ~ 10 23 cm" 2 ; the total 
molecular mass is a few 10 8 M Q . 

Accounting for high UV fluxes and variations in ki- 
netic temperature and assuming that the observed emis- 
sion arises from photon dominated regions (PDRs) re- 
solves the problems related to an LVG treatment of the 
radiative transfer. Spatial densities are as in the LVG case 
(n(H2) ~ 10 3 7 cm -3 and ~10 3 cm~ 3 for the high and low 
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excitation component, respectively), but 12 CO/ 13 CO in- 
tensity ratios <^10 indicate that the bulk of the CO emis- 
sion arises in UV-illuminated diffuse cloud fragments of 
small column density (A^H^) ~ 5 10 20 cm _2 /kms _1 ) and 
sub-parsec cloud sizes with area filling factors 3>1. Thus 
CO arises from quite a different gas component than the 
classical high density tracers (e.g. CS, HCN) that trace 
star formation rates more accurately. The dominance of 
such a diffuse molecular interclump medium also explains 
observed high [Ci]/CO line intensity ratios. PDR models 
do not allow a determination of the relative abundances of 
12 CO to 13 CO. Ignoring magnetic fields, the CO emitting 
gas appears to be close to the density limit for tidal dis- 
ruption. Neither changes in the 12 C/ 13 C abundance ratio 
nor variations of the incident far-UV flux provide good fits 
to the data for simulations of larger clouds. 

A warm diffuse ISM not only dominates the CO emis- 
sion in the starburst region of M82 but is also ubiq- 
uitous in the central region of our Galaxy, where tidal 
stress, cloud-cloud collisions, shocks, high gas pressure, 
and high stellar densities may all contribute to the for- 
mation of a highly fragmented molecular debris. 12 CO, 
12 CO/ 13 CO, and [C i]/CO line intensity ratios in NGC 253 
(and NGC 4945) suggest that the CO emission from the 
centers of these galaxies arises in a physical environment 
that is similar to that in M82. Starburst galaxies at 
large distances (z ~ 2.2-4.7) show 12 CO line intensity 
ratios that are consistent with those observed in M82. 
PDR models should be applicable to all these sources. 
12 CO/ 13 CO line intensity ratios 3>10, sometimes observed 
in nearby ultr aluminous mergers, require the presence of 
a particularly diffuse, extended molecular medium. Here 
[Cl]/CO abundance ratios should be as large or even 
larger than in M 82 and NGC 253. 
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sitions, only a few CO 6-5 spectra were presented from 
nearby galaxies (Harris et al. 1991; Wild et al. 1992). 

We have used the Heinrich-Hertz-Telescope (HHT) on 
Mt. Graham (Baars & Martin 1996) to map M82 in the 
CO J= 7-6, 4-3, and 13 CO 3-2 transitions. These data 
are complemented by new J = 2-1 and 1-0 spectra taken 
with the IRAM 30-m telescope. 



1. Introduction 

Low lying rotational transitions of CO are widely used 
as tracers of molecular hydrogen and are essential to de- 
termine dynamical properties and total molecular masses 
of galaxies. The widespread use of CO J = 1-0 and 2-1 
spectroscopy is however not sufficiently complemented by 
systematic surveys in higher rotational CO transitions to 
confine the excitation conditions of the dense interstellar 
medium (ISM). While the J = 1 and 2 states of CO are 
only 5.5 and 17 K above the ground level, the J = 3 to 7 
states are at 33, 55, 83, 116, and 155 K and trace a com- 
ponent of higher excitation. 'Critical densities', at which 
collisional deexcitation matches spontaneous decay in the 
optically thin limit, are ~ 10 5 ~ 6 cm~ 3 for CO J = 3-2 
to 7-6 in contrast to 10 35 and 10 43 cm -3 for the ground 
rotational CO transitions. 

Starburst galaxies are known to contain large amounts 
of molecular gas that may be heated to Tk; n ~ 100 K 
by young massive stars, cosmic rays or turbulent heat- 
ing. Therefore highly excited CO transitions, observed at 
submm- wavelengths, are the appropriate tool to study this 
interstellar gas component. Among the three nearest (D 
~ 3 Mpc) nuclear starburst galaxies, NGC 253, NGC 4945, 
and M 82 (NGC 3034) M 82 is most readily accessible from 
telescopes of the northern hemisphere. Containing one of 
the brightest IRAS point sources beyond the Magellanic 
Clouds (Sioo^im ~ 1000 Jy), M82 has been observed at a 
variety of wavelengths, ranging from the radio to the 7-ray 
domain of the electromagnetic spectrum. The starburst in 
M82 is likely triggered by a tidal interaction with M81, 
causing a high infrared luminosity (Lfir ~ 41O 1O L0), a 
high density of supernova remnants, and copious amounts 
of dense gas with strong OH and H2O masers and a 
large number of molecular high density tracers (for CO 
maps, see Sutton et al. 1983; Olofsson & Rydbeck 1984; 
Young & Scoville 1984; Nakai et al. 1986, 1987; Lo et al. 
1987; Loiseau et al. 1988, 1990; Phillips & Mampaso 1989; 
Turner et al. 1991; Tilanus et al. 1991; Sofue et al. 1992; 
White et al. 1994; Shen & Lo 1995; Kikumoto et al. 1998; 
Neininger et al. 1998). 

So far, few CO 4-3 maps of external galaxies were 
published (for M 51, M 82, M 83, and NGC6946 see White 
et al. 1994; Petitpas & Wilson 1998; Nieten et al. 1999). 
Among these M82 is the only true starburst galaxy but 
its CO 4-3 map (White et al. 1994) is confined to the very 
central region. With respect to higher rotational CO tran- 



2. Observations 

2.1. Observations with the Heinrich-Hertz-Telescope 

13 CO 3-2 (331GHz ~ 907/xm), 12 CO 4-3 (461 GHz ~ 
650/im), and 12 CO 7-6 (807 GHz ~ 372^m) line emis- 
sion was observed at the HHT during Feb. 1999 with 
beamwidths of ~ 22", 18", and 13", respectively. For 
the CO 3-2 and 4-3 transitions, SIS receivers were em- 
ployed; the CO 7-6 line was observed with a Hot Electron 
Bolometer (HEB) kindly provided by the Center for Astro- 
physics (Kawamura et al. 1999). The backends consisted 
of two acousto optical spectrometers, each with 2048 chan- 
nels (channel spacing ~ 480 kHz, frequency resolution ~ 
930 kHz) and a total bandwidth ~ 1 GHz. 

Spectra were taken using a wobbling (2 Hz) secondary 
mirror with a beam throw of ±120 to ±240" in azimuth. 
Scans obtained with reference positions on either azimuth 
were coadded to ensure fiat baselines. Receiver temper- 
atures were at the order of 170 K at 331 GHz, 150 K at 
461 GHz, and 1000 K at 807 GHz; system temperatures 
were ~ 900, 3500, and 8000 K on a Tj£ scale, respectively 

The receivers were sensitive to both sidebands. Any 
imbalance in the gains of the lower and upper sideband 
would thus lead to calibration errors. To account for this, 
we have observed Orion-KL (at 807 GHz) and IRC+10216 
(at 807 and 331 GHz) prior to our M82 measurements 
with the same receiver tuning setup. Peak temperatures 
were 70, 20, and 2.3 K on a T£ scale, respectively (cf. Howe 
et al. 1993; Groesbeck et al. 1994). At 461 GHz, Orion-KL 
(T£ ~ 70 K; cf. Schulz et al. 1995) was mapped but the 
tuning was later changed by AV = 200 km s -1 for M82 
(see Sect. 3.2). 

All results displayed are given on a main beam bright- 
ness temperature scale (T m b). This is related to T£ via 
T mb = T* A (F cS /B cS ) (see Downes 1989). Main beam ef- 
ficiencies, _B c ff, were 0.5, 0.38, and 0.36 at 330, 461, and 
806 GHz, as obtained from measurements of Saturn. For- 
ward hemisphere efficiencies are 0.9, 0.75, and 0.70, re- 
spectively (D. Muders, priv. comm.). With an rms surface 
accuracy of ~ 20^m (A/18 at 806 GHz), the HHT is quite 
accurate. This reduces the effect of the source coupling 
efficiency on the measured source size. 

At the center of M82, CO lineshapes depend sensi- 
tively on the position observed, so that in each of our maps 
the dynamical center could be identified with an accuracy 
better than 5". While relative pointing errors should be 
small when compared to the spacing (10") of our 13 CO J 
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Fig. 1. Intensities integrated over the velocity inter- 
val 25 to 385 km s -1 towards the nuclear starburst 
region of M82. 

Upper left: 13 CO J = 3-2 map over the central 50" x 
50". Contour levels are 14.2 to 71.0 by 7.1Kkms~ 1 . 
Upper right: 12 CO J — 4-3 map over the central 80" 
x 70". Contour levels are 124 to 620 by 62Kkms~ 1 . 
Lower left: 12 CO J = 7-6 map over the central 40" 
x 30". Contour levels are 64 to 320 by 32Kkms _1 . 
The nominal reference position is 0:1950 = 9 h 51 m 43 s , 
^1950 = 69° 55' 00", but in view of pointing errors 
(see Sects. 2.1 and 4.1), coordinates are best defined 
by a comparison with interfcromctric maps (e.g. Lo 
et al. 1987; Shen & Lo 1995; Neininger et al. 1998). 
For consistency with other studies, D = 3.25 Mpc is 
assumed throughout the paper. Recent distance de- 
terminations indicate slightly higher values (e.g. D = 
3.9±0.3Mpc; Sakai & Madore 1999). 



= 3-2 and CO 4-3 maps, relative pointing is less reliable 
in the case of the CO J = 7-6 map; here deviations may 
reach ±5" for a few positions. 



2.2. Observations with the IRAM 30-m telescope 

12 CO and 13 CO J = 1-0 and 2-1 observations were 
made with SIS receivers of high image sideband rejection 
(~25db for 1-0 and ^13 db for 2-1 line data) of the inner 
50" x 50" ( 12 CO; 5" spacing for the central 20", oth- 
erwise 10") and 20" x 20" ( 13 CO; 5" spacing) of M82 
in June 1999. The measurements were made in a posi- 
tion switching mode with the off-position displaced by 15' 
in right ascension. The two 12 CO lines (as well as those 
of 13 CO) were measured simultaneously. Beamwidths at 



115 (J = 1-0) and 230GHz (J = 2-1) were 21 and 13"; 
forward hemisphere and beam efficiencies were 0.92 and 
0.72 for the J = 1-0 and 0.89 and 0.45 for the J = 2-1 
data, respectively. Calibration was checked by observing 
IRC+10216. Measured line intensities were T m b = 17 and 
43 K for 12 CO J = 1-0 and 2-1 (channel spacings were 3.2 
and 1.6kms~ , respectively) and 2.3 and 6.0 K for 13 CO 
J = 1-0 and 2-1 (channel spacings: 2.7 and 1.4 km s -1 ; cf. 
Mauersberger et al. 1989). 



3. Results 



4 



R.Q. Mao, C. Henkel, A. Schulz et al.: The nuclear region of the starburst galaxy M82 



- C0(J=3-2) 6a - 50 : - 

: - Q 



IB 



Center 



(-10,-5); 




200 400 
Velocity (km/s) 



Fig. 2. Spectra towards the center (middle), the north- 
eastern (left), and the south-western (right) lobe. The 
temperature is given in units of main beam brightness 
temperature (K). Beamwidths are 13", 18", 22", 13", 13", 
21" and 21" for the 12 CO 7-6, 12 CO 4-3, 13 CO 3-2, 12 CO 
and 13 CO 2-1, and 12 CO and 13 CO 1-0 profiles, respec- 
tively. Offsets in arcsec relative to a nominal center po- 
sition (see Fig.|l|) are given in the upper right corner of 
each box. Small differences in these offsets are caused by 
pointing deviations (Sect. 2) or changes in source morphol- 
ogy (Sect. 4). All data are smoothed to a channel width of 
~10kms -1 . Linewidths (compare e.g. the J = 2-1 and 7-6 
profiles from the south-western lobe with the correspond- 
ing J = 3-2 and 1-0 spectra) are affected by differences 
in beamsizc. 
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Fig. 3. 13 CO J = 3-2 {upper panel) and 12 CO J = 4- 
3 (lower panel) channel maps. Center velocities are given 
in the upper right corner of each box. The radial veloc- 
ity range is 30kms _1 per image. The lowest black con- 
tour, the lowest white contour and the interval are re- 
spectively 1.4Kkms _1 , 7Kkms _1 , and 1.4Kkms _1 (up- 
per panel) and 14Kkms _1 , 70Kkms -1 , and 14Kkms _1 
(lower panel) on a T m b scale. 
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Fig. 4. Position-velocity diagram for CO J = 4-3 (solid 
lines; our data) and 2-1 (dashed lines; data from WeiB (in 
preparation)) along the line connecting the south-western 
and north-eastern lobe. Angular resolutions are 18" and 
13", respectively. The zero offset position is located half 
way between the emission peaks. CO contour levels are 
50%, 60%, 90%, and 99% of the respective peak inten- 
sities (see Table [l]) . 



3.1. Overall morphology of CO submillimeter line 
emission 

Maps of integrated 13 CO J = 3-2, 12 CO 4-3, and 12 CO 
7-6 intensity are presented in Fig.|l|. Detectable emission 
is strongly confined to the central part of the galaxy. In 
each map a particularly wide spectral feature, with slightly 
smaller peak intensity than the most intense lines in the 
SW and NE, could be identified with the dynamical cen- 
ter (see Fig.|| which also contains IRAM spectra). Two 
main peaks of emission are detected, being displaced by 
almost 10" from the center (this corresponds to a pro- 
jected distance of 150 pc); the south-western hot spot is 
most prominent, while evidence for a third peak in the NE 
(at ~ (17", 17")), seen in the CO J — 4-3 and 7-6 maps, 
is not conclusive. Velocity channel maps of the 13 CO J 
= 3-2 and 12 CO J = 4-3 line emission are displayed in 
Fig.||. These outline the extent of the emission at vari- 
ous velocities, show the dominant rotation pattern with 
the red-shifted lobe in the north-east and the blue-shifted 
lobe in the south-west, and indicate the rapid change in 
radial velocity near the dynamical center (cf. Neininger et 
al. 1998). 



3.2. Consistency of CO line temperatures 

Calibration at sub-millimeter wavelengths is critical be- 
cause of rapidly changing weather conditions, high atmo- 
spheric opacities, imbalances in the receiver gains between 
the sidebands, and because of uncertainties in beam and 
forward hemisphere efficiencies. Calibration uncertainties 
introduced by these effects may rise up to ±30% and a 
comparison with data published elsewhere is needed. In- 
tensities of our 13 CO J = 3-2 spectra from the center, 
the south-western, and north-eastern lobes are smaller 
by ~30% than those given by Tilanus et al. (1991; their 
Fig. 3). Our spectrum from IRC+10216 is however 30% 
stronger (on a T m b scale) than that given by Groesbeck 
et al. (1994), so that our scaling is intermediate between 
those of Tilanus et al. and Groesbeck et al. CO J = 4-3 
line intensities (Fig.|J) are ~30% larger than those given 
by White et al. (1994). Since their data were obtained 
with higher angular resolution (JCMT beamwidth: 11"), 
this difference is significant. A comparison with the three 
4-3 spectra from the 10-m CSO shown by Giisten et al. 
(1993) shows good agreement for the lobe positions. Dif- 
ferences in 4-3 peak temperatures between central and 
lobe positions are however less pronounced than reported 
by Giisten et al. (we find peak line temperature ratios of 
~1.3 instead of ~2.0 between the lobes and the center). 
Our CO 7-6 line intensities from Orion-KL agree well with 
those obtained by Howe et al. (1993) with the 10-m CSO 
antenna. 

Previous high angular resolution CO data from the 
30-m telescope (compare Fig. 9a of Loiseau et al. 1990 
with our Table were affected by uncertainties in the 
image sideband ratios. Our new data are not significantly 
affected by this problem, should be reliably calibrated 
within ±10%, and agree well with recent spectra obtained 
independently in an 'on-the-ffy' observing mode (WeiB, in 
preparation) . 

3.3. Line intensity ratios 

Claims that the large scale (<;20") integrated CO J = 2- 
1/1-0 line intensity ratio were much larger than unity (e.g. 
Knapp et al. 1980; Olofsson & Rydbeck 1984; Loiseau et 
al. 1990) can be firmly rejected. The more recently mea- 
sured ratios of 1.0 (Wild et al. 1992), 1.1 (Mauersberger 
et al. 1999) and 1.0-1.4 (Table 0) imply that 12 CO line 
intensities of the three lowest rotational CO transitions 
must be similar over the central 22". Our 13 CO J = 3-2 
line temperatures are a factor of ~10 smaller than those of 
the J = 3-2 12 CO transition observed with the same tele- 
scope (R. Wielebinski, priv. comm.). The situation with 
respect to the J = 6-5 line remains unclear. 6-5 spectra 
from the western lobe (Harris et al. 1991) and the center 
(Wild et al. 1992; their Fig. 2) do not allow a beam con- 
volution to the angular resolution of the lower frequency 
data and the beam pattern appears to be complex. Our J 
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Table 1. Line intensities (T = T m b), integrated line intensities (/ = jT nl ^dv), and line temperature ratios for 
a beamwidth of 22" a . CO J = 1-0 and 2-1 data were obtained with the IRAM 30-m telescope, higher excited 
rotational lines were measured with the HHT. 

Line or Line Ratio Position 5 

NE-Lobe Center SW-Lobe 



Line temperatures: 
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6.88±0.40 


7.40L0.36 






1.17L0.01 
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a) Given la errors in T m b and J T m bdf were determined from gaussian fits and neither include calibration uncertainties (see 
Sects. 2 and 3.2) nor errors caused by the limited spatial coverage of our CO 7-6 (Fig.0), 13 CO 2-1, and 13 CO 1-0 (Sect. 2) 
maps. For the beam convolution, gaussian distributions were assumed for the telescope beam and the source. While all 
measured positions were taken for beam convolution, the CO J = 4-3 beam size is almost as large as that for the 13 CO J = 
3-2 spectra. Therefore only five measured CO J — 4-3 spectra, the profile measured toward the specific position and those 
with offsets of 10" along the four cardinal directions, dominate the convolution. For the J = 7-6 line, spatially constraining 
the convolution to positions with weights >0.3 or A = 11 "7 (weight = exp[-4 In 2 (A 2 /(#bi - 9l 2 ))\, A: Offset; bl = 22", 6» b2 
= 18" or 13" for J = 4-3 and 7-6 spectra, respectively) modifies the result by only ~5%. 

b) For the offsets in arcsec, see FigM 

c) Because of the limited extent of the 13 CO J = 2-1 map (Sect. 2.2), 13 CO 2-1 values for the lobes should be considered with 
caution. The problem is less severe for the 13 CO J = 1-0 transition, because beam convolution is not required in this case. 



= 7-6 data show peak line temperatures of order 2-3 K. 
This is smaller than the 4 K measured in the J = 4-3 and 
lower J transitions with larger beam sizes. Beyond J = 
4-3 we thus find clear evidence for a weakening of CO 
emission with increasing rotational quantum number J. 



4. Spatial distributions: Are there differences? 

4.1. CO 

A comparison of the CO J = 4-3 data presented in Fig.[l] 
with those of White et al. (1994) shows a strong discrep- 
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ancy in the overall spatial distribution: Our CO J = 4-3 
map contains (at least) two maxima of emission, while the 
higher angular resolution data of White et al. (1994) have 
only one peak. Data with only one peak (that are based 
on spectra with sufficient resolution to separate the lobes) 
are also presented by Wild et al. (1992; their Fig. 10 and 
Table 5) for 13 CO and C 18 J =2-1. The original 13 CO 
J = 2-1 spectra displayed by Loiseau et al. (1988; their 
Fig. 1), however, clearly show a double-lobed distribution. 
Maps in the low-J 12 CO and 13 CO transitions (e.g. Lo 
et al. 1987; Nakai et al. 1987; Loiseau et al. 1988, 1990; 
Tilanus et al. 1991; Shen & Lo 1995; Neininger et al. 1998) 
as well as our J = 7-6 (Fig. 3), 2-1 and 1-0 data all show 
a double-lobed structure. We thus conclude that, in spite 
of previous evidence to the contrary, the overall spatial 
distribution of emission from highly excited CO shows two 
main centers of emission. 

Are the two main hotspots observed in the J = 1-0 and 
2-1 transitions identical with those observed in higher ex- 
cited CO lines? The position angles (east of north) of the 
lines connecting the hotspots are slightly smaller in our 
submillimeter data (Fig.^j) than in interferometric maps 
(cf. Lo et al. 1987; Shen & Lo 1995; Kikumoto et al. 
1987; Neininger et al. 1998). This is likely an error in our 
data caused by measurements at varying hour angles with 
pointing offsets along the azimuth and elevation axes. 

More significant is a difference in angular separation: 
While the two main peaks of line emission observed by us 
are separated by 15"±2", the interferometric maps (Lo et 
al. 1987; Shen & Lo 1995; Neininger et al. 1998) show a 
separation of 27"±2". This difference is larger than the 
positional uncertainties. A larger separation in the low-J 
transitions is also supported by filled-aperture measure- 
ments of Nakai et al. (1986, 1987) and Wild et al. (1992) 
for CO J = 1-0, by Loiseau et al. (1990) for CO J = 2-1, 
by Loiseau et al. (1988) for 13 CO J = 2-1, and by us for 
the 12 CO J = 1-0 and 2-1 lines. The CO J = 3-2 dis- 
tribution (Tilanus et al. 1991) shows an intermediate lobe 
separation. 

Which peaks are detected in the CO submillimeter 
lines? The two main lobes of CO emission are located al- 
most symmetrically with respect to the kinematical cen- 
ter, both in lines of low (e.g. Neininger et al. 1998) and 
high (see Fig. |J) excitation. We can therefore exclude that 
we see the north-eastern lobe and the 'compact central 
core' (see Shen & Lo 1995; Neininger et al. 1998) that are 
separated by ~15". Instead, the sudden drop of angular 
separation from 27" to 15" must reveal inhomogeneities 
in the molecular ring that are characterized by changes in 
density and temperature. In addition to the 'high CO exci- 
tation component' there exists a 'low CO excitation com- 
ponent', mainly emitting in the CO J = 1-0 and 2-1 lines. 
The transition in lobe separation occurs at the J = 3-2 
line: In CO J — 3-2 (Tilanus et al. 1991) the separation 
is still ~20". 13 CO J = 3-2 emission with smaller optical 



Table 2. Lobe separations ' 



Tracer 


Separation 


Reference 




arcsec 




Continuum 






2. 2 urn 


9±1 


McLeod et al. (1993) 6 


2.2/xm 


8±2 


Dietz et al. (1986) 


10.0/xm 


9±1 


Dietz et al. (1989) 


10.8Atm 


9±1 


Telesco et al. (1991) 


1 2 4/;m 


9±1 


Telesco Oezari (^<)<)2) 


19.2/wn 


11±2 


Telesco et al. (1991) 


19.5/im 


10±2 


Dietz et al. (1989) 


30.0/Mn 


10±2 


Telesco et al. (1991) 


450^m 


16±2 


Alton et al. (1999) 


450^m 


16±3 


Smith et al. (1991) 


450^m 


14±1 


Hughes et al. (1994) 


850/xm 


16±2 


Alton et al. (1999) 


1300^m 


18±4 


Kriigel et al. (1990) 


2700^m 


17±2 


Neininger et al. (1998) 


3400Aim 


18±4 


Seaquist et al. (1998) 


Spectral lines c 






Br 7 2.16/im 


10±2 


Lester et al. (1990) 


Br 7 2.16/im 


10±2 


Larkin et al. (1994) 


[Neil] 12.8/xm 


15±1 


A & L (1995) d 


CO 7-6 0.37 mm 


15±3 


This paper 


[Cl] 0.6 mm 


21±6 


White et al. (1994) e 


CO 4-3 0.65 mm 


14±3 


This paper 


CO 3-2 0.9 mm 


20±5 


Tilanus et al. (1991) 


13 CO 3-2 0.9 mm 


13±3 


This paper 


CO 2-1 1.3 mm 


24±3 


Loiseau et al. (1990) 


13 CO 2-1 1.4mm' 


15±5 


Loiseau et al. (1988) 


CO 1-0 2.6 mm 


26±2 


Shen & Lo (1995) 


CO 1-0 2.6 mm 


25±3 


Carlstrom (1988) 


CO 1-0 2.6 mm 


27±3 


Lo et al. (1987) 


CO 1-0 2.6 mm 


16±5 


Nakai et al. (1987) 


13 CO 1-0 2.7 mm 


27±2 


Neininger et al. (1998) 


CS 2-1 3.1mm 


26±3 


Baan et al. (1990) 


HCO+ 1-0 3.4 mm 9 


20±2 


Carlstrom (1988) 


HCN 1-0 3.4 mm 


24±3 


Carlstrom (1988) 


HI 21 cm 


29±3 


Weliachew et al. (1984) 



a) Errors are mainly by-eye-estimates. 

b) From Fig. 10 of McLeod et al. (1993) 

c) Lobe separations refer to maps displaying integrated line 
intensities. 

d) Achtermann & Lacey (1995) 

e) Since the CO J = 4-3 map of White et al. (1994) lacks the 
characteristic double lobed structure, their [C i] 492 GHz map 
should also be confirmed by independent measurements. 

f) The given lobe separation refers to Fig. 2 of Loiseau et al. 
(1988), displaying integrated line intensities. Peak line temper- 
atures (their Fig. 1) show a separation of ~25". 

g) The interferometric HCO + map of Carlstrom (1988; angu- 
lar resolution: 10") shows two pronounced peaks. It is however 
difficult to identify the lobes in the integrated intensity map of 
Seaquist et al. (1998; angular resolution 3'.'50 x 3'.'25). This 
may be caused by missing flux in the latter map. 
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depths and less photon trapping requires, however, higher 
excited gas so that the lobe separation becomes smaller. 

Since lobe separations are often comparable to the an- 
gular resolution of the particular observation, integrated 
intensity maps may be misleading as they easily exhibit 
structure dominated by the superposition of components 
with identical lines of sight but different velocities. To 
clarify the situation, we thus present in Fig.|| position- 
velocity maps of the CO J = 4-3 (our HHT map with 
highest signal-to-noise ratio and best relative pointing) 
and CO 2-1 (Weifi, in preparation) line emission. Beam 
widths are 18" and 13", respectively. Surprisingly, the p-v 
diagrams do not reproduce the strikingly different lobe sep- 
arations seen in the integrated intensity maps. This also 
holds when comparing CO J = 2-1 with 7-6 emission. 
This hints at wider line profiles for the higher excited CO 
transitions at the inner edges of the CO J = 1-0 and 2-1 
lobes. This is corroborated by Fig. [| that shows, for the 
CO J = 4-3 line, a slightly smaller lobe separation and 
a significantly larger full-width-to-half-power line width 
towards the inner edge of the south-western lobe. 

4-2. Other tracers of the interstellar medium 

Table |^ displays lobe separations determined by various 
tracers of atomic line, molecular line, or dust contin- 
uum emission. There are three preferred angular distances: 
~26", 15", and 10". Spectral lines of low excitation show 
two maxima at an angular distance of ~26". Highly ex- 
cited lines and the far-infrared and mm-wave continuum 
show a lobe separation of ~15". The continuum in the 
near and mid infrared and the Br 7 line show an angular 
spacing of ~10". 

Since the dynamical center of the galaxy is located 
not too far from the mid point of the line connecting the 
lobes, these data show evidence for an even more complex 
structure than indicated by the CO data alone. The hot 
dust and Br 7 emission from the inner parts of the ring 
may trace the most recent star formation activity as it 
propagates outwards into the molecular lobes located at 
larger galactocentric distances (Alton et al. 1999). 

Since the far infrared and submillimeter continuum are 
associated with the central portion of the ring, this is the 
location where column densities of the cool dense gas must 
be largest. 



5. Excitation analysis 

The J = 7-6 transition is the CO line with highest rota- 
tional quantum number so far observed in M 82. Measure- 
ments of this line widen the range of excitation accessible 
by CO data considerably. To analyse these excitation con- 
ditions, we have performed radiative transfer calculations 
using a Large Velocity Gradient (LVG) model describing 
a cloud of spherical geometry (see Appendix A) . 



5.1. Physical parameters: Our data 

The three submillimeter CO transitions measured by us 
arise from the high excitation component, show a similar 
spatial distribution, and can therefore be used simulta- 
neously in a radiative transfer analysis. To obtain true 
line intensity ratios, the CO 4-3 and 7-6 spectra were 
smoothed to the angular resolution of the 13 CO 3-2 data. 
Table [TJ displays these ratios for a beamwidth of 22", as- 
suming gaussian beamshapes. As indicated by the spatial 
distributions (see Fig.[l]), line ratios are similar toward the 
lobes and the central beam. CO 7-6/4-3 (integrated) line 
intensity ratios are at the order of 0.3, while 12 CO 4- 
3/ 13 CO 3-2 ratios are 7-10 (see Table 0). 



In Appendix A (Fig. A.l) we plot our LVG line in- 
tensity ratios as a function of density, kinetic tem- 
perature, 12 CO/ 13 CO ratio, and CO 'abundance' (i.e., 
A(CO)/gradF with Z(CO) = [CO]/[H 2 ]). For the line 
parameters given in Table [j], reasonable solutions can 
be found for small CO abundances (X(CO)/gradV^ ~ 
10- 6 pc/kms _1 ) and high 12 CO/ 13 CO ratios (~ 75). High 
CO abundances (X(CO)/gradV - 10~ 3 pc/kms" 1 ) and 
small 12 CO/ 13 CO ratios are however not providing a 
realistic solution, because H 2 densities would approach 
10 2 cm -3 and would thus become prohibitively small. Cal- 
culations assuming a plane-parallel instead of a spherical 
cloud geometry would yield even smaller densities. 

A comparison of lincwidths and velocity drifts with 
the size of the region suggests gradV^ ~ lkms _1 /pc. 
With this value and AT (CO) - 10" 5 - 10~ 4 (e.g. Blake 
et al. 1987; Farquhar et al. 1994) we can exclude a 
12 CO/ 13 CO abundance ratio as small as 25 and are 
guided to the four diagrams of Fig. |A.l| with log X (CO) 
~ log(A:(CO)/[gradT//kms" 1 pc- 1 ]) = -5 or -4 and 
12 CO/ 13 CO = 50 or 75. Our parameters are summarized 
in Table Ej. Kinetic temperatures are T^ n ~ 60 - 130 K, 



densities n(H 2 ) ~ 10 



3.3-3.9 



cm , and area filling factors 



/a,22" ~ 0.04 - 0.07 and f aA5 „ ~ 0.07 - 0.11. While 
the density is well constrained, the kinetic temperature 
is less well determined. Solutions with Tk; n ^ 150 K are 
however unlikely because they would require H 2 densi- 
ties <510 3 cm~ 3 . Even excluding such extreme solutions, 
the density of the CO emitting gas is small when com- 
pared with that of the Orion Hot Core (e.g. Schulz et 
al. 1995; van Dishoeck & Blake 1998), but column den- 
sities are large: With AT(CO) = 3.0810 18 [fico/cm" 3 ] 
[(AF/grady)/pc] cm~ 2 , we obtain the 22" beam averaged 
and cloud averaged column densities displayed in Table || 
For a line-of-sight source size of 350 pc and A^(H 2 ) 22 </ = 
10 23 cm -2 (see Table ||), the beam averaged mean molec- 
ular density is <n(H 2 )> 22 » ~ 200 cm~ 3 and the volume 
filling factor becomes / v ,22" = <^(H 2 )>/n(H 2 ) ~ 0.05. 
Since the ratio / v , 22" //a, 22" denotes the line of sight di- 
mension of the clouds in units of the beam size, we ob- 
tain with r cloud = 0.5tg22" D pc [/ v ,22»//a,22»] ~ 150 pc 
a characteristic cloud radius (that will be discussed and 
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Fig. 5. Peak line temperatures as a function of rotational 
quantum number J ( J — > J — 1) for a beam size of 22" 
toward the dynamical center of M82. 12 CO: Open square: 
Wild et al. (1992); stars: Mauersberger et al. (1999); open 
triangle: R. Wielebinski (priv. comm.); circles: this paper. 
13 CO (crosses): this paper. Results from radiative trans- 
fer calculations are given for 12 CO and 13 CO. Solid lines: 
X(CO)/gradF = 1CT 5 and 1.33 10" 7 pc/kms"\ n(H 2 ) 
= 8000cm" 3 , T kin = 70K; dashed lines: A(CO)/gradV 
= 10" 5 and 210" 7 pc/kms~ 1 , n(H 2 ) = 5000cm" 3 , 
T kin = 100K; dotted lines: A(CO)/gradI/ = 10" 4 
and l.SSlO^pc/kms" 1 , n(H 2 ) = 2000cm" 3 , T kin = 
80 K; dash-dotted lines: X(CO)/gradF = 10" 4 and 
210" 6 pc/kms" 1 , u(H 2 ) = 2000cm" 3 , T kin = 130K. 



revised in Sects. 6.2 and 6.3). The total molecular mass is 
Mmol,22» ~ 1-7 10 s M Q . 

5.2. Physical parameters: All CO data 

So far, we have only analysed the 'warm CO component' 
of M82, exclusively seen in 13 CO J = 3-2 and higher ex- 
cited rotational transitions. To combine these results with 
data from lower J rotational transitions and to further 
elucidate the physical state of the gas, Figs. || - ^ show (in- 
tegrated) line intensities as a function of quantum number 
J. Calibration errors are at the order of ±10% in the J = 
1-0 and 2-1 lines and ±20% in the higher excited lines. 

Figs. [5] and ^| show (integrated) line temperatures for 
a 22" beam toward the dynamical center of the galaxy. 
Fig.Q displays line temperatures for a 15" beam toward 
the lobes and demonstrates that CO excitation is similar 



Fig. 6. Integrated line intensities as a function of rota- 
tional quantum number J (J — * J — 1) for a beam size 
of 22" toward the dynamical center of M82. 12 CO: Open 
squares: Wild et al. (1992); open triangle: Glisten et al. 
(1993); stars: Mauersberger et al. (1999); circles: this pa- 
per. 13 CO (crosses): this paper. While most of our spa- 
tially convolved integrated line intensities were obtained 
from relatively extended maps (see Fig.|l|), the convolution 
of our 13 CO J = 2-1 data may be affected by insufficient 
spatial coverage (see Sect. 2.2). The CO J = 2-1, 3-2, and 
4-3 data from Wild et al. (1992) and Gusten et al. (1993) 
are based on few observed positions only. These cover the 
ridge of strong CO emission and only allow to estimate 
intensities for the dynamical center of M82. To account 
for the missing positions, the resulting integrated intensi- 
ties were multiplied by 0.85 (the factor was deduced from 
the CO J = 3-2 map of Tilanus et al. 1991). Results from 
radiative transfer calculations are also given (for the pa- 
rameters, see Fig.jHI). 



toward the south-western and north-eastern hotspot. Re- 
sults from radiative transfer calculations (see Sect. 5.1 and 
Appendix A) are also given. Input parameters correspond 
to the four boxes in Fig. A.l that provide a promising fit 

= —5 or 



to the data (log(A(CO)/[grady/kms 1 pc~ 
-4 and 12 CO/ 13 CO = 50 or 75). 

Apparently there is a problem with the CO J = 4- 
3 line: The integrated intensity (Fig.^J) does not allow a 
reasonable fit, while the peak intensity (Fig. ||) is 'appro- 
priate'. This is caused by the narrow lineshape of our 4-3 
spectrum (Fig.|^). Compared with other lines, the inte- 
grated line intensity is too small, but the peak line temper- 
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Fig. 7. Peak line temperatures as a function of rotational 
quantum number J (J — * J — 1) for a beam size of 
15" toward the south-western (squares) and north-eastern 
(stars) lobe of M82. All temperatures were rescaled for 
coupling the 15" beam to a source solid angle (at a given 
velocity) of Q = lOOarcsec 2 (e.g. Lord et al. 1996). 12 CO 
J = 1-0: Nakai et al. (1987; 16"), this paper (22"); J = 

2- 1: Wild et al. (1992; 13"), this paper (13"); J = 3-2: 
Tilanus et al. (1991; 14"); J = 4-3: Giisten et al. (1993; 
15"), this paper (18"); J = 7-6: this paper (13"). 13 CO 
(crosses): J = 1-0 and 2-1: this paper (22" and 13"); J = 

3- 2: Tilanus et al. (1991; 14"), this paper (22"). Results 
from radiative transfer calculations are also given (for the 
parameters, see Fig.||). 



ature is almost 'normal'. The CO J = 4-3 profile, shown 
by Giisten et al. (1993) for the central position, is also 
weak, both with respect to its peak and integrated line 
intensity. The emission from the lobes, however, fits into 
the general trend (Fig.|?]). An interpretation in terms of a 
diminished lobe separation for CO J > 4-3 (see Sect. 4.1) 
is not conclusive: The 22" beam centered on the dynam- 
ical core of M 82 is mainly confined to the inner parts of 
the lobes and should not be greatly affected by emission 
from further out. A CO J = 4-3 deficiency is neither seen 
in Fig.]?] nor in a corresponding plot showing integrated 
intensities for a 15" beam. 

To summarize: Data from the J = 4-3 line are contra- 
dictory so that convincing evidence for a true anomaly 
is missing. In spite of differences in lobe separation 
(Sect. 4.1), the CO data can be reproduced, within ob- 
servational errors, with densities, temperatures, filling 



nuclear region of the starburst galaxy M 82 
Table 3. LVG model parameters (compare with Table ||[) 



Footnote 

High excitation comp. 



12 CO/ 13 CO 


50-75 


a 


Tkin 


60-130 K 




n(H 2 ) 


10 3.3...3.9 cm -3 




fa, 22" 


0.04 ... 0.07 


b 


/a, 15" 


0.07 ... 0.11 


b 


f 

Jv,22" 




c 


r(CO J = 1-0) 


0.5 - 4.5 




r(CO J = 2-1) 


3.5 - 15 




r(CO J = 4-3) 


10 - 40 




r(CO J = 7-6) 


6-25 




N(CO) 22 „ 


-5 10 18 cm" 2 


d 


iv(H 2 )22" 


~10 23 cm" 2 


d 


JV(CO)doud 


~10 20 cm- 2 


e 


Af(H2) c loud 


~10 24 - 25 cm" 2 


e 




~10 22 cm- 2 /kms _1 




<n(H 2 )>22» 


-200 cm -3 




^cloud 


~150 pc 


f 


A^mol,22" 


~ 1-7 10 s M 




Low excitation comp. 






n(H 2 ) 


~10 3 cm -3 





a) Abundance ratio 

b) / ai22 " : Area filling factor for the central 22" 

c) / Vj 22" : Volume filling factor for the central 22" 

d) iV(CO)22", -A(H2) 2 2": 22" beam averaged total column 
densities 

e) -/V(CO) c i ou d, A r (H2)cioud: 22" cloud averaged total column 
densities 

f) 7" c i ud: Characteristic molecular cloud size 



factors, [CO]/[H 2 ] abundances, and 12 CO/ 13 CO isotope 
ratios deduced from the three submillimeter transitions 
mapped by us with the HHT (Fig.|l|). 

There are few constraints for the low excitation compo- 
nent that is mainly seen in the CO J = 1-0 and 2-1 lines. 
The column density must be smaller than for the high 
excitation component because the latter is more closely 
related to the far infrared and submillimeter continuum 
from the dust (Sect. 4.2). In view of the remarkable num- 
ber of 'super'-star clusters (O'Connell et al. 1995) and 
supernovae (Kronberg et al. 1985) near the outer portions 
of the ring, cloud temperatures for the low excitation com- 
ponent should also be T kin > 10 K. For T kin - 50 K, LVG 
densities are at the order of n(H 2 ) ~ 10 3 cm~ 3 or less. 

6. Discussion 

6.1. A comparison with other LVG simulations 

The most detailed models of CO emission from M 82 were 
so far provided by Giisten et al. (1993). In order to fit the 
12 CO/ 13 CO J = 1-0 and 2-1 line intensity ratios then 



11 



available, they rejected a one component LVG scenario 
and introduced two gas components, one of low (n(H2) ~ 
10 3 cm -3 ) and one of high (~10 5 cm -3 ) density. The low 
density component is similar to that proposed by us for 
the gas mainly emitting in the CO J — 1-0 and 2-1 lines 
(see Sect. 4.1). Our parameters for the high CO excitation 
component agree to a large extent with those of the one 
component scenario of Giisten et al. (1993). 

For the highly excited gas we thus also find opti- 
cally thick 12 CO low- J emission and a high [ 12 CO]/[ 13 CO] 
abundance ratio (^50; see Table ||), that is further sup- 
ported by an independent chain of arguments involving 
CN and 13 CN data (see Henkel et al. 1998). While the LVG 
model result is shown to be inconclusive in Sect. 6.3.2, the 
CN data suggest that the [ 12 CO]/[ 13 CO] ratio is larger 
than that observed in the galactic center region. Likely 
explanations are radial infall of 13 CO deficient gas from 
the outer parts of the galaxy or a 12 C excess in the ejecta 
from massive stars (e.g. Henkel & Mauersberger 1993). 
In agreement with Giisten et al. (1993) we also find that 
excitation in the south-western and north-eastern lobe is 
similar. Since our new data allow us to constrain kinetic 
temperatures to the high end of those predicted by the one 
component scenario of Giisten et al. (Tki n <; 30-70 K), our 
densities are at the low end (n(H2) Ss 10 4 cm -3 , since the 
product Tkin ^.(H^) 1 / 2 is approximately conserved among 
models simulating optically thick subthermal CO emis- 
sion). There is a remarkable agreement with respect to 
CO column density, H2 density, and molecular gas mass 
between various studies (cf. Tilanus et al. 1991; Wild et 
al. 1992; Giisten et al. 1993). 

6.2. Intrinsic inconsistencies of the model 

Our CO data could be reproduced by assuming the pres- 
ence of two gas components. Selecting such a 'best' model, 
we could discriminate between previously proposed sce- 
narios and could constrain cloud conditions giving rise to 
highly excited CO emission to slightly higher tempera- 
tures and smaller densities than previously suggested. On 
a deeper level, however, there remain problems. While cal- 
culated CO column densities (Sect. 5.1) appear to be cor- 
rect (see e.g. Fig. 2 of Lo et al. 1987 and Fig. 3 of Smith 
et al. 1991), an obvious puzzle is the large volume filling 
factor (/ v ,22" ~ 0.05) that is comparable to the area filling 
factor (Sect. 5.1). This forced us to postulate a character- 
istic cloud radius (r c i ou d ~ 150 pc) that encompasses a 
large part of the studied volume. Such a large cloud ra- 
dius is inconsistent with the expectation of similar scale 

lengths along the three dimensions (/ v ,22" ~ /a 22" ) ana - 
with the spatial fine structure deduced from CO (e.g. Shen 
& Lo 1995; Neininger et al. 1998), high density tracers like 
HCN and HCO+ (e.g. Brouillet & Schilke 1993; Paglione 
et al. 1997; Seaquist et al. 1998), and infrared fine struc- 
ture lines (e.g. Lugten et al. 1986; Lord et al. 1996). 



Another problem is the density of the gas observed. A 
density of n(H2) 10 4 cm~ 3 is small when compared to 
the prototypical 'hot core' associated with the Orion neb- 
ula, violating a sometimes noted similarity between these 
tiny, highly obscured galactic star forming regions and the 
more extended starbursts studied in external galaxies (e.g. 
Lo et al. 1987; Wolfire et al. 1990). More seriously, the 
molecular density determined by us for the starburst in 
M 82 is smaller than most theoretical studies and obser- 
vational data permit: Assuming 'reasonable' density stuc- 
ture (n(H2) ex r -1 "" ~ 2 ) and accounting for the intense UV 
field, Brouillet & Schilke (1993) find that molecular clouds 
with densities less than a few times 10 4 cm~ 3 should not 
exist in the central region of M 82. For the transition re- 
gion between the atomic and molecular gas they propose 
a density of 10 4 ~ 5 cm~ 3 which disagrees with the range of 
densities deduced from our LVG analysis in Sect. 5.1 (for 
cloud stability against tidal stress, see Appendix B). 

Studies of high density tracers (e.g. CS, HCN, and 
HCO + ) commonly reveal densities n(H2) ^ 10 4 cm -3 (e.g. 
Mauersberger & Henkel 1989; Brouillet & Schilke 1993; 
Paglione et al. 1997; Seaquist et al. 1998), even for the low 
excitation component (compare Mauersberger & Henkel 
1989 with Baan et al. 1990). In view of T dust = 48 K 
(Hughes et al. 1994; Colbert et al. 1999), radiative excita- 
tion by infrared photons should not significantly alter the 
density estimates (see Carrol & Goldsmith 1981 and Ap- 
pendix C). Millimeter wave recombination lines indicate 
the presence of an ionized component with low filling fac- 
tor and high electron density (n c > 10 45 cm~ 3 ; Seaquist 
et al. 1996). This gas may be associated with a popula- 
tion of (ultra)compact HII regions or with shock ionized 
dense molecular material. Far infrared fine-structure lines, 
if tracing the interface between the molecular and the ion- 
ized gas, indicate densities of 10 3 3-4 cm -3 (e.g. Lugten 
et al. 1986; Wolfire et al. 1990; Lord et al. 1996; Stutzki 
et al. 1997; Colbert et al. 1999) that may agree too well 
with those derived by us for CO (Sect. 5.1). 

An interesting aspect is also provided by kinetic tem- 
perature estimates: With Tki n ~ 20-60 K (Seaquist et al. 
1998; this temperature is consistent with the estimated 
cosmic ray flux; see Volk et al. 1989, Suchkov et al. 1993, 
and Fig. 2 of Farquhar et al. 1994) the temperature of the 
dense molecular gas appears to agree fairly well with that 
of the dust (T dust = 48 K; Hughes et al. 1994; Colbert 
1999). This is further supported by an apparent lack of 
CH3OH and SiO emission, two tracers of high temperature 
gas that are easily seen in other nearby starburst galaxies 
(Henkel et al. 1991; Mauersberger & Henkel 1993). The 
surrounding neutral and ionized layers have larger tem- 
peratures, at the order of 50-100 to 200 K (e.g. Lugten et 
al. 1986; Wolfire et al. 1990; Lord et al. 1996; Colbert et al. 
1999). An analysis of [Ci] 492 and 809 GHz emission from 
the south-western lobe of M82 (Stutzki et al. 1997) shows 
a particularly striking similarity in density, temperature, 
and area (but not volume) filling factor with our LVG 
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parameters derived for CO (Sect. 5.1; the [Ci]/CO abun- 
dance ratio then becomes ~0.4). It seems as if CO were 
an integral part of the dense atomic gas layers of M82. 

6.3. A PDR model for the CO emission from M82 
6.3.1. General aspects 

To resolve inconsistencies related to CO fine scale struc- 
ture, H.2 density, and kinetic temperature (Sect. 6.2), we 
note that a significant fraction of the molecular gas in the 
Milky Way lies in photon dominated regions (PDRs; see 
Hollenbach & Tielens 1997). A PDR model is also suc- 
cessfully applied to simultaneously explain CO and [Cn] 
line intensities in the central region of the late-type spiral 
IC 342 that is not a starburst galaxy but that is believed 
to be a face-on 'mirror image' of the galactic center region 
(Schulz et al. in preparation). 

In the central region of M82 with its high UV flux, 
the bulk of the interstellar line radiation likely arises from 
PDRs. As we have seen (Sect. 6.2), the coolest gas compo- 
nent is observed toward the dense cloud cores, a situation 
that is consistent with the inside-out temperature gradi- 
ent expected in the case of PDRs. While CO excitation by 
infrared radiation from the dust can be neglected (see Ap- 
pendix C), the impact of UV photons may strongly affect 
rotational level populations and measured CO line inten- 
sities in regions with high UV flux. Estimating the far-UV 
flux (6.0 - 13.6 eV) in the starburst region, Wolfire et al. 
(1990), Stacey et al. (1991), Lord et al. (1996), and Colbert 
et al. (1999) find x ~ lO 28-39 incident far-UV flux 
in units of the local galactic flux, 1.6 10 -3 ergs -1 cm -2 ). 
Such a high value leads to outer cloud layers that are 
predominantly heated by collisions with electrons photoe- 
jected from dust grains or through collisional deexcitation 
of vibrationally excited, UV-pumped H 2 . 

Exploring effects of finite cloud size for clumps with 
plane-parallel geometry, Koster et al. (1994) presented 
comprehensive computations of CO rotational line inten- 
sities from PDRs. In contrast to the constant temperature 
LVG treatment presented in Sect. 5, PDR models account 
for strong kinetic temperature gradients. The low- J lines 
typically arise from deeper inside the clouds than the mid- 
J lines observed by us (Sect. 3.1). Reproducing observed 
line ratios of optically thick 12 CO emission, PDR model 
densities are larger than those derived from a one temper- 
ature, one density LVG model. PDR line intensity ratios 
only surpass unity if the higher- J line has a critical density 
(for the values, see Sect. 1) that is smaller than the actual 
gas density. Then, both lines are approximately thermal- 
ized but the higher J-line is emitted from warmer layers 
further out. For densities at the order of 5 10 3 cm -3 (the 
density deduced from our one component LVG model in 
Sect. 5.1), 12 CO PDR line temperatures rapidly decrease 
with rotational quantum number J. For \ ~ 10 3 and 
N(H 2 ) ~ 10 22 cm -2 /kms -1 (Sect. 5.1), densities are then 



at the order of 10 4 and 10 5 cm -3 for the low and high 
excitation components, respectively. 

For the gas with low CO excitation, a density of n.(H2) 
~ 10 4 cm -3 is sufficiently high to explain the detection of 
CS J = 2-1 emission (e.g. Mauersberger & Henkel 1989; 
Baan et al. 1990). A density of n(H 2 ) ~ 10 5 cm -3 in the 
high excitation region would fulfill all theoretical (Brouil- 
let & Schilke 1993) and observational density requirements 
outlined in Sect. 6.2. Furthermore, the volume filling factor 
/v,22" = <n(H2)>22"/ n (H2) would drop by one to two or- 
ders of magnitude below the value estimated in Sect. 5.1. If 
the area filling factor is not drastically altered, this yields 
reasonable molecular cloud radii SjlO pc (see Sect. 5.1 and 
Lugten et al. 1986; Wolfire et al. 1990; Brouillet & Schilke 
1993; Shen & Lo 1995; Seaquist et al. 1996; Stutzki et al. 
1997). The apparently intermediate temperature of the 
CO emitting gas between those of the dust and the dense 
atomic medium is naturally explained by mid- J CO emis- 
sion predominantly arising in the heated surface layers of 
dense molecular clouds. 

To summarize: PDR simulations of 12 CO emission re- 
move inconsistencies related to spatial fine scale structure, 
density, and kinetic temperature. So far published PDR 
results fail however when 13 CO is also considered. Calcu- 
lated 12 CO/ 13 CO line intensity ratios are much smaller 
than the ratios observed (see Table [j] and Koster et al. 
1994) and a more detailed numerical analysis is therefore 
needed. 

6.3.2. PDR model calculations 

For a numerical approach we note that the beam averaged 
column density (^(112)22" ~ 10 23 cm -2 ) and the beam 
averaged density (<n(H 2 )>22" ~ 200cm -3 ) are observa- 
tionally determined (Sects. 5.1 and 6.2) and do not de- 
pend on the choice of the excitation model if a significant 
fraction of the dust is associated with molecular clouds. 
Complementing these boundary conditions, we obtain (cf. 
Sect. 5.1) 

r c ioud ,„ 200cm" 3 /n(H 2 )pDR 

[PCJ JCO,obscrvcd/-7cO,PDR 

and 

iy(H 2 )pD R 10 23 /300 

[cm- 2 kms- 1 ] TcO,observed/2cO,PDR 

with 200 cm -3 /n(H2)pDR denoting the volume filling fac- 
tor, Tco,observed/7co,PDR being the area filling factor, 
"(H2)pdr giving the average density, and A^H^pdr 
representing the average column density of an indi- 
vidual clump per kms -1 (300 km s -1 is the total CO 
linewidth of the nuclear region of M82; A ? (H2)pdr 
= (4/3) • r c i ou d • n(H 2 )pDR and n(H 2 )pDR-cioud-surfacc = 
0.5-n(H2)pDR for the assumed cloud geometry and den- 
sity structure, see below). 
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Since embedded clumps of dense gas may more often 
be spherical than plane-parallel, spherical clouds with a 
power law density distribution (n(r) ~ r" 1 - 5 ) were mod- 
eled with the dust being heated by the external UV ra- 
diation and intrinsic infrared emission (cf. Storzer et al. 
1996; 2000). While we can reproduce the observed relative 
intensities of the various rotational CO transitions with a 
density at the order of n(H.2) ~ 10 5 cm -3 (see Fig.^ and 
Sect. 6.3.1), r c i ou d and iV(H 2 ) are inconsistent with Eqs. 1 
and 2. Furthermore, computed 12 CO/ 13 CO line intensity 
ratios show little dependence on cloud structure and re- 
main with characteristic values of 2-4 (see also Gierens et 
al. 1992 for 12 C/ 13 C = 40; Koster et al. 1994 and Storzer et 
al. 2000 for 12 C/ 13 C = 67) much smaller than the observed 
ratios of 10-15 (see Table Varying \ over the permit- 
ted range (io 3 - 3±0 - 4 ) does not change these line intensity 
ratios significantly. Drastically increasing the carbon iso- 
tope ratio leads to 12 C/ 13 C > 100 which is inconsistent 
with data from the solar system, the Milky Way, and the 
Magellanic Clouds (e.g. Wilson & Rood 1994; Chin et al. 
1999). 

Other parameters that can be varied are the den- 
sity and column density of the far-UV irradiated cloud. 
The 12 CO lines are formed in warm layers near the sur- 
face whereas 13 CO lines are predominantly emitted from 
cooler regions deeper inside. With the column density be- 
ing fixed, a higher density leads to a moderate increase in 
the 12 CO/ 13 CO line intensity ratio (see Fig.||), because 
12 CO is then emitted from CO layers closer to the heated 
cloud surface. To reproduce observed 12 CO/ 13 CO line ra- 
tios in this way requires, however, densities ^10 8 cm -3 . 

More important are variations in column density: A 
spherical UV illuminated cloud is effectively smaller in 
13 CO than in 12 CO and this difference becomes more pro- 
nounced when iV(Ha) decreases (see Storzer et al. 2000). 
In the extreme case of totally photodissociated 13 CO, 
some 12 CO may still exist in the clump so that, in prin- 
ciple, 12 CO/ 13 CO line intensity ratios up to very large 
values can be reproduced. Once the 13 CO abundance gets 
very small, 12 CO also tends to become optically thin. In 
Fig.^j we show expected line temperatures in this opti- 
cally thin regime, with n(H-2) = 5 10 3 cm -3 successfully 
reproducing measured 12 CO line intensity ratios for the 
high CO excitation component. 12 CO/ 13 CO = 19, 10, 8.5, 
9.4, 13, 20, and 30 for J = 1-0 ... 7-6, respectively. This 
is compatible with observed ratios of 15, 10, and 9 for 
the three lowest rotational transitions. For the low excita- 
tion component, n(H2) ~ 10 3 cm -3 ; Fig. ^| shows a fit with 
12 CO/ 13 CO = 17, 10, and 11 for J = 1-0 to 3-2. 

It is remarkable that PDR model densities are simi- 
lar to those derived with our LVG simulation, contradict- 
ing theoretical predictions of minimum densities in excess 
of 10 4 cm" 3 (Brouillet & Schilke 1993). Column densities 
(iV(H2) ~ 5 10 20 cm^/kms" 1 ) and line temperatures of 
individual cloudlets are, however, drastically different. An 
average cloudlet shows optically thin, not optically thick 



Table 4. PDR model parameters 



Footnote 



High excitation comp. 






12 CO/ 13 CO 


undetermined 


a 


n(Ha) 


~10 3 ' 7 cm -3 




/a, 22" 


~ 15-200 


b 


/v,22" 


~ 0.05 


c 


iV(CO) 22 » 


-510 18 cm -2 


d 


iV(H2)22» 


~10 23 cm* 2 


d 


A r (H 2 )cloud 


~5 10 20 cm^/kms" 1 


c 


<n(H 2 )>22" 


-200 cm" 3 




^cloud 


~0.025pc 


f 


A^ mo l,22" 


~ 1-7 10 s M 




Low excitation comp. 






n(H 2 ) 


~10 3 cm" 3 




/a, 22" 


~15 


b 


^"cloud 


~0.15pc 


f 


,) 12 CO/ 13 CO: abundance ratio 



b) / a ,22" : Area filling factor for the central 22"; the lower 
value is derived from a comparison of modelled and observed 
line intensities, the higher value results from a comparison of 
total 22" beam averaged and individual cloud column density. 

c ) /v,22" : Volume filling factor for the central 22" (see 
Sect. 5.1) 

d) iV(CO) 2 2", iV(H2) 2 2": 22 " beam averaged total column 
density 

e) iV(H2)cioud: Column density of an individual cloud 

f) fcioud: Characteristic molecular cloud size 



CO line emission. Line temperatures of individual clumps 
are much smaller than the observed average over the in- 
ner 300 pc, so that the PDR area filling factor / a is not 
■Cl as in the LVG approximation but 3>1. Thus not a 
higher density (as anticipated in Sect. 6.3.1) but a higher 
/ a leads to plausible r c i ou( j values in the sub-parsec range 
(see Eq. 1). Both / a 3> 1 and r cloud < lpc are consis- 
tent with early CO studies, favoring CO lines of low op- 
tical depth. These conclusions are supported by observa- 
tions of fine-structure lines in the infrared (e.g. Olofsson 
& Rydbeck 1984; Lugten et al. 1986; Wolfire et al. 1990; 
Brouillet & Schilke 1993; Schilke et al. 1993; Lord et al. 
1996). In view of our PDR simulations, the high density 
(>10 cm -3 ) ionized gas found by Seaquist et al. (1996) 
is likely related to (ultra)compact HII regions and not to 
evaporating dense molecular clouds. 

Obviously, our study does not imply that the core of 
M 82 does not contain regions of high density and col- 
umn density. The bulk of the CO emission, however, must 
arise from a warm low density interclump medium. The 
gas may be barely dense enough to avoid tidal disruption 
(see Appendix B). For typical cloud cores in the galactic 
disk, Eq. 5 of Larson (1981) predicts a column density of 
7V(H 2 ) - 10 22 cm" 2 . With iV(H 2 ) ~ 5 10 20 cm^/kms" 1 
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(Table |J) this infers a linewidth of 20kms _1 . The rela- 
tions defined for approximately virialized clouds by Lar- 
son (1981) then yield cloud sizes well in excess of 100 pc, 
in strong contradiction with our PDR model. The bulk of 
the CO emission therefore arises from gas that may not 
be virialized. We conclude that CO is tracing a different 
gas component than molecular high density tracers like CS 
or HCN. The dominance of CO emission from a diffuse 
medium with cloud fragments of low column density per 
kms -1 , X = N(H.2)/Igo conversion ratio, and CO inten- 
sity explains why HCN is a better tracer of star formation 
and infrared luminosity (Solomon et al. 1992a). It also 
explains why [Cl]/CO abundance ratios are higher than 
those observed in spatially confined star forming regions 
of the galactic disk (see also Sect. 6.2, Schilke et al. 1993, 
and White et al. 1994). 

An argument against our scenario could be that 
Eqs. 1 and 2 only hold approximately: The denominator 
Tco,obscrvod/Tco,PDR is too small by an order of magni- 
tude to yield the proper r c i ou( j value assumed in the PDR 
approximation, while it is too large by the same amount 
to provide the proper -/V(H2)pdr value. In order to check 
whether our scenario is plausible, we therefore extend our 
analysis to other galaxies. 

6.4- Other galaxies 

6.4.1. The Milky Way and other 'quiescent' galaxies 

Studying the physics of the molecular gas in the galac- 
tic center region, Dahmen et al. (1998) compared 12 CO 
and C 18 J = 1-0 data obtained with a linear resolu- 
tion of ~22pc (8'). Lineshapes, spatial distributions, and 
line ratios indicate the presence of an extended diffuse 
gas component that is not apparent in studies focusing on 
individual cloud cores (see also Oka et al. 1998a). A signif- 
icant fraction of the gas is not virialized. Higher resolution 
data (e.g. Huttemeister et al. 1998) reveal that the dense 
molecular gas (n(H2) iS 10 4 cm -3 ) is relatively cool (Tki n 
~ 25 K) in comparison to that at lower density (<~ 100 K 
at a few 10 3 cm~ 3 ). Apparently, a warm diffuse molecu- 
lar medium is not only ubiquitous in the starburst galaxy 
M82 but also dominates the CO emission in the more 
quiescent central region of the Milky Way (for extended 
CO maps of a galactic disk star forming region, see Wil- 
son et al. 1999). This also holds for the nuclear regions of 
IC 342 (Schulz et al. in preparation) and NGC 7331 (Israel 
& Baas 1999), another two galaxies with Lfir ~ 10 10 L Q . 
Like M82, NGC 7331 shows a warmer inner and cooler 
outer dust ring and a molecular ring that appears to be 
related to its cool dust component. 

Shocks and cloud-cloud collisions induced by the pres- 
ence of bars (e.g. Achtermann & Lacey 1995; Morris & 
Serabyn 1996; Fux 1997, 1999; Huttemeister et al. 1998), 
tidal disruption of clouds near the center (Glisten 1989), a 
high gas pressure (e.g. Heifer & Blitz 1997) that can help 
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Fig. 8. PDR model peak line temperatures as a func- 
tion of rotational quantum number J (J — > J - 1) for 
clouds with high column density. Dotted lines with open 
squares: n(Ha) = 210 4 cm~ 3 , the curve with high Tt, dis- 
plays 12 CO, the lower one 13 CO; dashed lines with as- 
terisks: n(H2) = 2 10 5 cm -3 ; solid lines with open circles: 
n(H 2 ) = 2 10 6 cm" 3 . For all models, N(R 2 ) = 2 10 22 cm" 2 , 
implying cloud radii of 7.5 10 17 , 7.5 10 16 , and 7.5 10 15 cm, 
respectively. Assumed linewidth of an individual cloud and 
its carbon isotope ratio: 1.2 kms~ and 67. For the defini- 
tion of n(H2) = n(H2)pDRj see Sect. 6.3.2. 



to keep the gas molecular, and a high stellar density that 
can affect molecular cloud dynamics (e.g. Mauersberger et 
al. 1996b; Oka et al. 1998b) may all contribute to the dis- 
integration of molecular clouds and to the synthesis of an 
extended warm molecular spray consisting of low column 
density cloud fragments. 

6.4.2. Nearby starburst galaxies 

So far, M 82 is the only starburst galaxy where we could 
show that a CO-LVG excitation analysis (inferring f v ~ 
/ a ; see Sect. 5.1) does not lead to results that are free of 
inconsistencies. The lack of CO J = 7-6 data does not 
constrain LVG parameters sufficiently to search for a sim- 
ilar discrepancy in the other two nearby starburst galaxies 
NGC 253 and NGC 4945. This implies that previous stud- 
ies did not provide sufficient motivation to replace LVG 
excitation analyses by PDR scenarios. 

Towards NGC 253, M82, and NGC 4945 12 CO and 
12 CO/ 13 CO line intensity ratios are similar. This even 
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Fig. 9. PDR model peak line temperatures as a func- 
tion of rotational quantum number J (J — ► J — 1) 
for low column density clouds, reproducing line inten- 
sity ratios for the 'high' and 'low' CO excitation com- 
ponents. Dotted lines with open squares ('low exci- 
tation'): 7i(H 2 ) = 10 3 cm -3 , r c i ou d = 4.510 17 cm, and 
7V(H2) = 6 1 20 cm -2 ; dashed lines with asterisks ('high 
excitation'): n(H2) = 510 3 cm -3 , r c i ouc i = 7.510 16 cm, 
and A^(H2) = 5 10 20 cm -2 ; solid lines with open circles: 
n(H 2 ) = 210 4 cm" 3 , r c i oud = 1.510 16 cm, and N(R 2 ) 
-1 10 20 cm -2 . Assumed linewidth of an individual cloud and 
its carbon isotope ratio: 1.2 km s -1 and 67, respectively. 
For the definition of n(Ha) = n(H-2)pT>K> see Sect. 6.3.2. 



holds for absolute integrated line intensities (within 30%) 
for a 22-24" beam size (compare Table [j] with Mauers- 
berger et al. 1996a, b; Harrison et al. 1999). Not sur- 
prisingly, published LVG densities and area filling fac- 
tors for NGC253 and NGC4945 match those for M82 
(Henkel et al. 1994; Mauersberger et al. 1996a, b; Harrison 
et al. 1999). Following the procedure outlined in Sect. 5.1, 
average densities along the line-of-sight (<n(H2)> ~ 
100 cm~ 3 ) yield volume filling factors at the order of / v ,22" 
= <rt(H2)>/?i(H2) ~ 0.02. Since volume filling factors 
seem to be slightly smaller than area filling factors, av- 
erage cloud sizes become ^30 pc (for the Milky Way, see 
Oka et al. 1998b), not enough to request the use of PDR 
models. 

Because of the similarity of CO line strengths and line 
ratios and since high [Ci]/CO ratios are observed toward 
both M 82 and NGC 253 (Schilke et al. 1993; White et al. 
1994; Harrison et al. 1995; Israel et al. 1995; Stutzki et al. 



1997) we feel nevertheless that PDR simulations relating 
the bulk of the CO emission to a warm diffuse molecular 
medium are relevant not only to M 82 but to NGC 253 and 
NGC 4945 as well. 

6.4.3. Mergers 

Active disk galaxies exhibit 12 CO/ 13 CO ratios ~ 10 (this 
is independent of inclination, Hubble Type, and metallic- 
ity; Sage & Isbell 1991). Perturbed gas-rich mergers with 
infrared luminosities <; 10 11 L Q tend to show larger val- 
ues, up to 40-60 (Henkel & Mauersberger 1993). This may 
be caused by excitation effects (see Aalto et al. 1999) or 
by a deficiency of 13 CO emission since Lfir is better cor- 
related with 12 CO (Taniguchi & Ohyama 1998). Inflow of 
gas from the outer disks may provide 13 CO deficient gas 
and nucleosynthesis in shortlived massive stars may fur- 
ther enhance 12 C (Casoli et al. 1992). While this may help 
to enhance 12 CO/ 13 CO line intensity ratios, this is not suf- 
ficient to explain measured 12 CO/ 13 CO values within the 
context of PDRs. Instead we suggest that the direct in- 
teraction of two galactic nuclei and their associated disks 
is even more efficient than the presence of a bar to trig- 
ger cloud-cloud collisions and to create a warm diffuse 
molecular debris containing a large number of small cloud 
fragments. If such fragments are smaller than in individual 
non-merging galaxies, higher 12 CO/ 13 CO line ratios might 
result. While the detailed spatial CO fine-structure may 
only be revealed by next generation mm-wave telescopes 
(e.g. Downes 1999), this would imply that in mergers with 
12 CO/ 13 CO > 10, [Ci]/CO line intensity ratios should be 
as large as or even larger than in M 82 and NGC 253. For 
a first detection of [C i] in a merger, see Gerin & Philips 
(1998). 

6.4.4. Galaxies at high redshifts 

Is it a general property of the integrated spectrum of a 
starburst that, beginning with the CO J = 4-3 transi- 
tion, line intensities (accounting for beam dilution) grad- 
ually decrease with increasing rotational quantum number 
J? Ground based measurements of nearby galaxies in mid 
and high-J CO transitions require exceptional weather 
conditions. In galaxies of high redshift, however, many 
otherwise inaccessible transitions are shifted into the ob- 
servable mm- and submm- wavelength bands. While lin- 
ear resolutions remain poor, the bulk of the CO emis- 
sion from higher excited states should arise from the nu- 
clear starburst environment. Suitable examples of distant 
sources with a number of detected CO transitions are 
IRAS F 1024+4724 (z = 2.286), the Cloverleaf quasar (z 
= 2.558), and BR 1202-0725 (z = 4.692). Toward IRAS 
F 1024+4724, the J = 3-2 line is stronger than the 4-3 and 
6-5 lines; line intensity ratios w.r.t. 3-2 are ^0.75 and 0.6 
(Solomon et al. 1992b). Toward the Cloverleaf quasar, CO 
lines are characterized by a constant or increasing bright- 
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ness temperature from 3-2 to 4-3, followed by constant 
or (more likely) decreasing line intensities in the higher J 
transitions (Barvainis ct al. 1997). Toward BR 1202-0725, 
the J = 7-6 to 5-4 line intensity ratio is ^0.65 (Omont 
et al. 1996). Although the sample is too small for a re- 
liable statistical analysis, it seems that CO line intensity 
ratios are fairly uniform in starbursts and depend little on 
redshift (age) and temperature (T cm b/[K] = 2.73 (1+z)) 
of the microwave background. Unfortunately, 12 CO/ 13 CO 
line intensity ratios are not yet known. If they are <T0, 
densities at the order of a few 10 3 cm -3 and small column 
densities (^10 21 cm~ 2 /kms -1 ) should be a characteristic 
feature for the bulk of the CO emitting cloudlets in all 
starburst galaxies. 

7. Conclusions 

We have studied millimeter and submillimeter CO line 
emission up to the J = 7-6 rotational transition toward 
the central region of the starburst galaxy M 82 and obtain 
the following main results: 

(1) The spatial structure of the millimeter and submillime- 
ter CO emission is distinct. While integrated intensity 
maps suggest that the lobe separation of the low- J 
transitions is ^26", it is ~15" for the mid-J transi- 
tions. Major-axis position-velocity maps in the CO J 
= 2-1 and 4-3 lines show however agreement in the 
lobe positions. This indicates that, at the inner edges 
of the low-J CO lobes, line profiles are wider in the 
higher excited CO transitions. We thus distinguish be- 
tween a 'low' and a 'high' CO excitation component, 
the latter coinciding with the main source of millimeter 
and submillimeter dust emission. 

(2) An LVG excitation analysis of CO submillimeter 
lines leads to internal inconsistencies. While mea- 
sured line intensities are reproduced with T^ n ~ 60 
- 130K, n (Ha) ~ 10 3 - 3 ~ 3 ' 9 cm -3 , cloud averaged col- 
umn densities iV(CO) c i ou d ~ 10 20 and A r (H 2 ) c ioud ~ 
10 24-25 cm -2 ; pC0]/[ 13 C0] abundance ratios ^50, 
and a total molecular mass of a few 10 8 M Q , area filling 
factors (/ a ~ 0.05-0.10) and volume filling factors (/ v 
~ 0.05) are similar. This results in cloud sizes that do 
not match their angular scale. On the other hand, the 
resulting H2 column density is consistent with that de- 
rived from the dust continuum at millimeter and sub- 
millimeter wavelengths. For the low excitation compo- 
nent, densities are ~ 10 3 cm~ 3 . 

(3) An application of PDR models resolves the inconsis- 
tencies of the LVG calculations. LVG densities, column 
densities, and total mass are confirmed. The bulk of 
the CO emission arises, however, from a diffuse, low 
column density (iV(Ha) ~ 5 10 20 cm _2 /kms _1 ) inter- 
clump medium with small X — A^H^/Ico conver- 
sion factors, area filling factors 3>1, and sub-parsec 
cloud sizes. The relations denned by Larson (1981) are 



not fulfilled and the gas may not be virialized. Such 
a scenario explains why CS or HCN are better trac- 
ers of global star formation rate and infrared lumi- 
nosity than CO. Our scenario also explains observed 
high [C i]/CO line intensity ratios, while relative abun- 
dances of 12 CO versus 13 CO cannot be accurately de- 
termined. Higher column density clouds, even account- 
ing for variations in far-UV flux and 12 C/ 13 C isotope 
ratios, do not reproduce observed 12 CO/ 13 CO line in- 
tensity ratios ^10. Densities are close to the minimum 
values required for tidal stability in the absence of mag- 
netic fields. 

(4) In regard to 12 CO line intensity ratios, the central re- 
gion of M 82 appears to be representative for the en- 
tire family of starburst galaxies, both at small and at 
high rcdshifts. A comparison of the starburst regions 
in M 82 and NGC 253 demonstrates that this similarity 
extends to 12 CO/ 13 CO and [Ci]/CO line intensity ra- 
tios. The large 12 CO/ 13 CO line intensity ratios (>10) 
observed toward 'nearby' mergers prove, however, that 
differences exist at least w.r.t. rare CO isotopomers. 
Galaxy pairs with such high 12 CO/ 13 CO line ratios 
require the presence of a particularly diffuse highly 
fragmented low column density ISM. 

Apparently, dropping the assumption of constant tem- 
perature in the CO excitation model is a necessary step to 
provide a self-consistent approach to the physical proper- 
ties of molecular clouds in the nuclear starburst region of 
M 82. While the use of PDR models is crucial for a better 
understanding of the molecular gas phase in a starburst 
environment, important information is still missing. In- 
terferometic observations of high density tracers (e.g. CN, 
CS, HCN, HNC, N 2 H+), coupled with PDR model calcu- 
lations including chemical aspects, are needed to fully un- 
derstand the spatial morphology, density distribution, and 
molecular excitation of this archetypical starburst com- 
plex. An interesting aspect is provided by HCO + J — 1-0 
line emission (Table ||). The bulk of this emission might 
arise from regions intermediate between those of the low 
and high CO excitation component. Since this molecule 
(as well as N 2 H + ) is a sensitive tracer of ionization condi- 
tions in the dense gas, a detailed knowledge of its spatial 
distribution would be crucial for a better understanding 
of structure and excitation. 

So far, models were calculated for the cosmic ray flux 
of the solar neighbourhood. A flux enhancement by two 
to three orders of magnitude (with all necessary chemical 
implications) has still to be incorporated into PDR codes 
(but see Schilke et al. 1993). Another important quantity 
is the spatial distribution of the UV flux. For M82, we 
do not know the variation of the UV flux as a function of 
galactocentric radius. 
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Appendix A: The radiative transfer model 

Systematic motions and microturbulence are frequently 
used as simplifying assumptions to facilitate treatment of 
line formation in molecular clouds. For the large amplitude 
systematic motions assumed in the Large Velocity Gradi- 
ent (LVG) approximation, the source functions are locally 
defined. For the homogeneous velocity field assumed in 
the microturbulent approximation, the source functions 
are generally coupled throughout the cloud by scattered 
radiation. The crucial difference between the two approxi- 
mations is that between local and non-local excitation and 
both can be viewed as limiting cases for the treatment of 
molecular line formation. In the case of CO, line intensi- 
ties determined with a standard LVG model (e.g. Castor 
1970; Scoville & Solomon 1974) do not differ by more than 
a factor of three from models using the microturbulent 
approach (White 1977; see also Ossenkopf 1997). This is 
within the uncertainties that can be attributed to cloud 
geometry. Furthermore, Wild et al. (1992) found no signif- 
icant difference between results from their 'clumpy cloud' 
and standard LVG models. Given the inhomogeneity of 
the ISM, the assumption of uniform physical conditions is 
crude and only yields average gas properties. In view of 
the limited quality of the available data and the lack of 
information on source morphology and its fine-scale struc- 
ture, however, an LVG treatment of the radiative transfer 
is an appropriate first step to analyse M82. 

Applying our LVG model to simulate a cloud of spheri- 
cal geometry, CO collisional cross sections were taken from 
Green & Chapman (1978). For Tk; n > 100 K, we assumed 

1 /2 

collision rates Cy oc T kin . Using instead the collision rates 
recommended by De Jong et al. (1975; these include higher 
temperatures) leads to similar results. A comparison of 
J=2-l/J=l-0 line intensity ratios with those of Castets 
et al. (1990), that were computed with yet another set of 
collision rates, also shows consistency within 20%. 



Fig. |A1| displays our calculated CO 7-6/4-3 and CO 
4-3/ 13 CO 3-2 line intensity ratios as a function of density 
(10 2 - 7 cm- 3 ), kinetic temperature (5-150K), 12 CO/ 13 CO 
abundance ratio (25-75), and CO fractional abundance 
in terms of A(CO)/grad V (10~ 6 -- 3 pc/kms" 1 ; A(CO) 
is the fractional abundance parameter, gradV denotes 
the velocity gradient in kms _1 /pc). The choice of 
A(CO)/gradV is motivated by a source size of ~ 300 pc 
and a velocity range of ^300 km s^ 1 (this leads to gradU 
~ lkms~ pc -1 ), a solar system [C]/[H] abundance ra- 
tio of 3.5 10~ 4 (Grevesse et al. 1994), and the assumption 
that a significant fraction of the available carbon (^10%) 
is forming CO molecules. Since molecular clouds with sig- 
nificant CO 7-6 emission must be warm, strong CO frac- 
tionation in favor of enhanced 13 C abundances (cf. Watson 



disk; e.g. Henkel et al. 1985), and 12 CO/ 13 CO = 75 (close 
to the carbon isotope ratio of the local ISM and the solar 
system; e.g. Stahl et al. 1989; Stahl & Wilson 1992). 

Appendix B: Cloud stability 

In order for a cloud to be gravitationally bound, it must 
be sufficiently dense to withstand the tidal stresses caused 
by the gravitational potential of the galaxy. Neglecting 
rotation, turbulence, and magnetic fields, and compensat- 
ing tidal forces by a cloud's own gravity, we can derive 
a minimum density that is required for survival in a hos- 
tile medium. With V m t ~ 140 km s -1 at a galactocentric 
radius of R il 75 pc from the dynamical center of M82 
(Neininger et al. 1998) and applying Eq. 5 of Giisten and 
Downes (1980), we find a limiting cloud density of n m i n 
~ 10 3 - 2 - 3 - 5 [i?/120pc]- 2 cm- 3 . R = 120pc refers to the 
molecular lobes of M82, located at offsets ±7'.'5 from the 
dynamical center of the galaxy (see Table |J). The mini- 
mum density n m ; n is near the lower limit of the densities 
deduced in Sect. 5.1 with an LVG model and just matches 
densities derived in Sect. 6.3.2 with a PDR model. 

Appendix C: Radiative pumping 

Because the kinetic temperature is constrained by the CO 
J=7-6/J=4-3 line temperature ratios (Table [l]), colli- 
sional excitation to excited vibrational or electronic levels, 
^ 3000 K above the ground state, is not effective. Densi- 



ties and temperatures obtained with Fig. A.l might how- 



et al. 1976) can be excluded. Fig. A.l thus presents cal- 
culations for 12 CO/ 13 CO = 25 (the carbon isotope ratio 
in the galactic center region; e.g. Wilson & Rood 1994), 
12 CO/ 13 CO = 50 (the 12 C/ 13 C ratio in the inner galactic 



ever be affected by radiative excitation, that is neglected 
in our LVG treatment (apart of the 2.7 K microwave back- 
ground) . 

To operate efficiently, radiative pumping must be 
faster than the collisional rates for rotational excitation 
in the v = vibrational ground state. With <av> ~ 
210~ 11 cm 3 s~ 1 (Green & Chapman 1978) and a density 
of 5 10 3 cm~ 3 , collision rates are at the order of 10 _7 s _1 . 
Direct rotational excitation by A ~ 1 mm photons from 
the dust is inefficient, because Tdust •C 1. For vibrational 
exitation by 4.7/im photons we obtain with the Einstein 
coefficients A mt = 6.1 10~ 6 and 3.4 lO^s" 1 for the J 
= 4-3 and 7-6 rotational transitions, A v= i^q ~ 20 s -1 
(Kirby-Docken & Liu 1978), and the condition Tdust ~ 
3Q7QK/ln[A v=1 ^oMrot] (see Eq.6 of Carroll & Gold- 
smith 1981; i/v=i->o ~ 6.410 13 Hz) a minimum dust tem- 
perature of Tdust ~ 200 K. This is much larger than Td us t 
= 48 K, estimated by Hughes et al. (1994) and Colbert 
et al. (1999). The condition for efficient pumping by in- 
frared photons is therefore, if at all, only fulfilled for a 
small part of the molecular complex in the central region 
of M82 (cf. McLeod et al. 1993). Adopting in spite of 
this Tdust ~ 200 K, the radiative pumping rate becomes 
J5 v= o^i M v =i-o ~ 10~ 15 s _1 (or less, if beam dilution plays 
a role) . This is negligible in comparison with the collisional 
pumping rate. 



18 



R.Q. Mao, C. Henkel, A. Schulz et al.: The nuclear region of the starburst galaxy M82 




2 4 6 



log n(H2) 



Fig. A.l. Large Velocity Gradient (LVG) model results for X(CO)/grady = 1(T 6 , 1(T 5 , 10~ 4 , and 10~ 3 pc/kms -1 
(X(CO) is the fractional CO abundance parameter) and 12 CO/ 13 CO abundance ratios of 25, 50, and 75. The logarithm 
of X(CO)/gradF and the 12 CO/ 13 CO abundance ratio are given for each box in the upper right hand corner. Solid 
lines: CO 7-6/4-3 line intensity ratios (contour levels from lower left to upper right: 0.01, 0.03, 0.1, 0.3, 0.5, 0.7, 0.8, 
and 0.9). Dotted lines: 12 CO 4-3/ 13 CO 3-2 line intensity ratios (contour levels from lower right to upper left: 1., 5., 
10., 15., and 25.; for 12 CO/ 13 CO = 25, line intensity ratios are <15, for 12 CO/ 13 CO — 50, line intensity ratios are 
<25). 
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